Dynamical Superfluid-Insulator Transition 
in a Chain of Weakly Coupled Bose-Einstein Condensates 
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We predict a dynammical classical superfluid-insulator transition (CSIT) in a Bose-Einstein con- 
densate (BEC) trapped in an optical and a magnetic potential. In the tight-binding limit, this 
system realizes an array of weakly-coupled condensates driven by an external harmonic field. For 
small displacements of the parabolic trap about the equilibrium position, the BEC center of mass 
oscillates with the relative phases of neighbouring condensates locked at the same (oscillating) value. 
For large displacements, the BEC remains localized on the side of the harmonic trap. This is caused 
by a randomization of the relative phases, while the coherence of each individual condensate in the 
array is preserved. The CSIT is attributed to a discrete modulational instability, occurring when 
the BEC center of mass velocity is larger than a critical value, proportional to the tunneling rate 
between adjacent sites. 
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The recent experimental investigations of the dynam- 
ical properties of a Bose-Einstein condensate trapped in 
optical potentials WflK have led to a rapidly growing 
interest in this topic . The spatial and temporal 

coherence of matter waves emitted at different heights 
of the gravitational field has been proven in after 
loading a condensate in a vertical optical trap. Num- 
ber squeezed (non-classical) states have been realized in 
1^ . In 1^ Bloch oscillations and interband transitions in 
an accelerating lattice have been observed. In the 
optical potential was superimposed on a harmonic mag- 
netic trap, realizing a chain of weakly coupled conden- 
sates (i.e., a Josephson junction array) driven by an ex- 
ternal parabolic field. For small initial displacements, the 
condensate center of mass oscillated symmetrically with 
the relative phases among adjacent sites locked togheter. 

Here we demonstrate that for large displacements, 
when the velocity of the center of mass reaches a crit- 
ical value (proportional to the tunneling rate beteen ad- 
jacent sites), the BEC abruptly stops on the side of the 
harmonic trap (i.e., without reaching its center). We de- 
fine an order parameter for the system, and show that 
this dynamical transition from a "superfluid" to an "in- 
sulator" regime is associated with a randomization of 
the relative phases among different wells. As wc will 
discuss below, this transition has a classical (mean-field 
or Gross-Pitaevskii) nature, and it is different (but with 
some analogies) from the quantum superfluid-insulator 
(Mott) transition caused by the number squeezing of 
the quantum states in each well p. It also differs from 
the Landau dissipation mechanism, occurring in (quasi- 
)homogeneous systems when the velocity of the conden- 
sate is larger than the sound speed Q|. Rather, the CSIT 
is driven by a modulational instability (MI) that causes 



an exponential growth of small perturbations of a carrier 
wave, as a result of the interplay between dispersion and 
nonlinearity. The MI is a general feature of discrete as 
well as continuum nonlinear wave equations. Its demon- 
strations span a diverse set of disciplines ranging from 
fluid dynamics (where it is usually referred to as the 
Benjamin-Feir instability) and nonlinear optics [ p2| to 
plasma physics [ p^ . One of the early contexts in which 
its significance was appreciated was the linear stability 
analysis of deep water waves. It was only much later rec- 
ognized that the conditions for MI would be significantly 
modified for discrete settings relevant to, for instance, the 
local denaturation of DNA iQ or coupled arrays of opti- 
cal waveguides In the latter case, the relevant model 
is the discrete nonlinear Scrodinger equation (DNLS), 
and its MI conditions were discussed in |l^]. In this let- 
ter we propose an experiment to observe a superfluid- 
insulator mean-field dynamical transition (as a conse- 
quence of the MI), with weakly coupled Bose-Einstein 
condensates driven by an external harmonic field. 

The BEC dynamics is governed by the Gross-Pitaevskii 
equation (GPE): 
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-, a is the s-wave scattering length, m 



where 

the atomic mass and /x the chemical potential. The con- 
densate wave function <&(r, t) is normalized to the total 
number of condensate atoms Nt^ and we consider a re- 
pulsive interatomic interaction a > 0. The external po- 
tential Vext is given by the sum of the harmonic confining 
potential Vjvf = ^[^^xX^- 
tential Vl = Vb cos^ {kx) 



'Uj'^{y'^ + z'^)] and the optical po- 



The valleys of the potential are 
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separated by a "lattice spacing" of A/2, with A = 27r/fc. 
We consider a chemical potential fj, << Vq and the trans- 
verse degrees of freedom to be frozen by a tight magnetic 
confinement, so as to justify the study of the system in 
an effective one-dimensional geometry. 

In the tight-binding approximation ^{r, t) = 
^JNt 4'nit)4'nir), with the wavefunction of the 
condensate in the j-th site of the array, weakly coupled 
in the barrier region with the wavefunctions 4'j±i of the 
condensates in the neighbour sites. It is then possible to 
map the GPE onto the DNLS §: 
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with K ~ ~ J rf?^[^(V0„ • V0„+i) + 0„Vext0n+i] pro- 
portional to the microscopic tunneling rate between adja- 
cent sites, U = goNr J dr(j)\ and e„ = / (if[|^(V(/)„)2 + 

(1)^ The DNLS Hamil- 
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with iip'^^'ipn canonically conjugate variables. Both Ji 
and the norm | ip^ p= 1 are integrals of the motion. 

Let us consider, first, the case e„ — (which corre- 
sponds to neglecting the effect of the harmonic trap). 
Stationary solutions of Eq. (^ are plane waves ip — 
'0oexp[i(fc7i — i/t)], of frequency = — 2_ftrcos(fc) + t/|V'oP- 
The stability analysis of such states can be carried out 
by perturbing the carrier wave with small amplitude 
phonons: tpn = V'oe'^'""''*) + we^(«"""*' + ^;*e~'(9"+'^*). 
The DNLS excitation spectrum (for e„ = 0) is then given 
by: 

uj± = 2i4rsin(A:) sin((7)± 
cos2(fc) sin''(|) + 2X[/|VoP cos(fc) sin2(|). (4) 

The carrier wave becomes modulationally unstable when 
the eigenfrequency uj in Eq.(^ becomes imaginary: 



C/ItAoP > -2i<rcos(fc)sin2(g/2), 



(5) 



namely, when cos(fc) < 0. Therefore, if the interatomic 
interaction is repulsive {U > 0), the system suffers an 
exponential growth of perturbations when 7r/2 < fc < 
37r/4. This result will remain valid in the case of non- 
homogeneous travelling wave-packets driven by external 
fields, when their width is much larger than the wave 
length associated with the collective motion. This con- 
clusion can be further understood in the light of the col- 
lective coordinate equations of motion developed in [Q. 
Generally speaking, the mapping of the GPE into the 
DNLS allows for the study of solitons and localized exci- 
tations as well as dynamical instabilities in the framework 



of the lattice theory |T7-19|. We remark, however, that 
the MI is a general feature of the GPE with a periodic 
external potential, and not necesserly in the tight bind- 
ing limit. In the perturbative limit, ^ » Vq, and in 
absence of external drivingficlds (i.e., with Vi\i ~ 0), the 
MI has been studied in [g|j8| (of course, in this limit the 
MI condition doffers from Eq.(^)). 

The effect of the exponential growth of phonon modes 
of arbitrary momenta in the DNLS leads to an effective 
dephasing among different sites of the lattice. Indeed the 
phases of each condensate enter into a "running regime" , 
with an angular velocity different from site to site and 
proportional to the local (on-site) effective chemical po- 
tential. The complete delocalization in momentum space 
leads to a strong localization in real space, hence to the 
appearance of localized structures of large amplitude (see 
Fig. |l|) . This localization has also been attributed (in the 
absence of any external potential e„ = 0) to the presence 
of the so-called Peierls-Nabarro barrier which pins 
such large amplitude solutions pl| , not allowing them to 
propagate. The excess kinetic energy is partially stored 
to high-frequency internal "ac" oscillations among ad- 
jacent wells (see also [^,|3|), and partially converted 
to wakes of small amplitude extended wave radiation 

The CSIT can be observed experimentally by condens- 
ing, firstly, the atomic gas in both the magnetic and the 
optical trap, and, then, adiabatically displacing the mag- 
netic field from its initial position. For small displace- 
ments, in line with the findings of the system coher- 
ently oscillates about the center of the potential. If we 
rewrite ipj = y^^"^' ' ^^^^ implies that the phase differ- 
ence between sites is given by (/)j+i(i) — (j)j{t) = Acj){t). 



The center of mass ^ 
A(j) will then satisfy 



'^^jnj and the phase difference 



h—iit) = 2K sinA</>(t) 
dt 

?i^A<^(i) = -2 17 e(t). 
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Eqs. (ph have the usual form of the Josephson equa- 
tions [P4| and indicate that the overall array of bosonic 
Josephson junctions behaves as a single Josephson junc- 
tion, whose critical current is iK/Ti. The collective co- 
herence was experimentally demonstrated in 0| by the 
interference pattern obtained upon releasing the conden- 
sates from the optical and magnetic traps. 

To monitor the dynamical transition of interest, we 
define < k >= X^fc^lV'feP = A0, with ipk the Fourier 
transform of the condensate wave-function. In the co- 
herent, small amplitude oscillations regime, the quasi- 
momentum < k > exhibits regular oscillations (see Fig. 
^). However, for < k >> 7r/2, the system becomes mod- 
ulationally unstable and localization ensues. The critical 
initial displacement ^cr can therefore be obtained from 
Eqs. (^) with Atp = Tr/2. In lattice units: 
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In Fig. ||we plot < fc > vs. time for three initial displace- 
ments. When £^{t — 0) is smaller than ^cr, the average 
momentum < k > oscillates in accordance with Eq. 
Note that for ^{t = 80) the system approaches very close 
to the instability line. When the initial displacement is 
larger than the critical value, < k > abruptly drops as 
soon as it crosses the critical point. This is accompa- 
nied by the sudden arresting of the BEG center of mass 
(cf. Fig. 1^) and by a collective dephasing (cf. Fig. ||). 
The key experimental signature would be the disappear- 
ance of the interference fringes after the expansion of the 
BEG, while the center of mass of the resulting cloud will 
be resting on the side of the trap's center. 

The values of the parameters in DNLS are Vq = 5Eji 

,2 

with A = 795 nm and the recoil energy Er = ^^^^ , 
Nt = 50000, K = 5.5- IQ-^Er and 17 = 1.5 • 10~^Er: 
the critical displacement from Eq. (0) is ^cr ~ 84 sites, 
in good agreement with our numerical findings (cf. Figs 
2 and ^). The loss of coherence is highlighted in Fig. 
3, where we plot the temporal evolution of the modulus 
squared of an effective complex order parameter measur- 
ing the overall coherence of the system, defined as: 



4' = E^.-V';+i- 



(8) 



When the collective oscillations are coherent, the value 
of the order parameter is — 1 (see the cases 

^{t = 0) = 40,80 in Fig. ^. On the other hand, 
a complete dephasing is characterized by |5'(t)P = 0, 
and occurs for ^(t ~ 0) > S,cr, or, equivalently, when 
< k >> 7r/2 (cf. Fig. 0). It is worth noting that such 
randomization takes place between the phases of conden- 
sates localized in different wells, even though each one 
remains internally coherent. This is shown in Fig. |^, 
where we plot the density for different times below (a,b) 
and above (c,d) the critical displacement. The numer- 
ical solutions of the DNLSE and the GPE are in good 
agreement. The motion of the center of mass in the su- 
percritical case is reported in Fig. ^ where the numer- 
ical solutions of the DNLS and the full one-dimensional 
GPE (|l|) are compared. In both cases the system stops 
at < X >~ 35 fim (with a slight difference between the 
DNLS and the GPE predictions) , while the center of the 
harmonic trap is located at x = fj,m. From Eqs. ^ 
we can calculate the critical current, i.e., the maximum 
allowed velocity in the coherent transmission of matter 
waves: by setting Ac/) = tt/2 we readily see that the crit- 
ical velocity of the center of mass, S^max is equal to the 
critical current per particle Ic = 2K/h. In dimensional 
units: 



DNLS numerical result, and close to the numerical GPE 
value Vc = 1.18 iim/ms. 

From the above findings, we can conclude that the ef- 
fect of the MI is to dephase the system. In the effec- 
tive ID geometry we have considered, such dephasing is 
strong enough to stop the falling condensate. In higher 
dimensions, the dephasing can be partial, and will only 
damp the BEG motion. Yet, its onset will still be given 
by Eq. (^ . The GSIT regards a classical field (the solu- 
tion of the GPE) and it is qualitatively different from the 
quantum Mott insulator-superfluid transition in meso- 
scopic Josephson junction chains, which is driven by the 
competition between zero-point quantum phase fluctua- 
tions and the Josephson coupling energy. Yet, it is pos- 
sible to draw an analogy. In the former GSIT case, the 
insulator regime is associated with a vanishing temporal 
correlation among the phases of each site, each phase still 
being meaningful in the GPE sense. The quantum transi- 
tion is also driven by a loss of phase correlations induced 
by the localization of atoms in each site, which, however, 
arises from the non-commuting nature of the number- 
phase observables. Glearly, such quantum fluctuations 
cannot be captured within the GPE framework. Also, 
the latter transition is reversible (i.e., long-range phase 
coherence is restored upon adiabatically increase of the 
tunnel coupling), while the former is not. It is worth 
noting that very recent experimental works [ p6[ have il- 
lustrated the existence (and reversibility) of the quantum 
phase transition, rendering the experimental verification 
of the classical dynamic transition suggested herein, a 
natural next step for experimental studies. In conclu- 
sion, we notice that the MI can be studied in term of one 
(or several) bifurcation points in an effective stationary 
Hamiltonian. Such bifurcation points separate regions 
with different symmetries, and it is common in the liter- 
ature to study such dynamical transitions in terms of an 
order parameter ||2^ , borrowing the language and con- 
cepts of statistical phase transitions. This mapping, in 
the specific case of MI, deserves further studies. 

The modulational instability (and the consequent 
superfluid-insulator transition) studied here can also be 
observed with different experimental setups (e.g., with 
the condensate at the center of the harmonic trap and 
with the laser moving across). In fact, similar MI and 
pinning results have been obtained in the case in which 
the harmonic trap is displaced at a constant speed ex- 
ceeding a critical value These results illustrate the 
generality and importance of the effects of the MI mech- 
anism in the motion of Bose-Einstein condensates and 
underscore its potential in inducing localization and de- 
phasing of such coherent structures. 
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FIG. 1. The density calculated at different times 0, 20, 
40 ms (from the right of each figure to the left) as a func- 
tion of the position, with initial displacements ^(0) = 50, 120 
sites, which are, respectively, below and above the critical 
value xi^r « 84 (0). The GPE (b,d), and the DNLS (a,c) 
wavefunctions normalized to 1 are compared. (a,b) ^(0) = 50 
sites; (c,d) ^(0) = 120 sites. The external parabolic potential, 
which drives the oscillations, is centered at a; = 0. 
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FIG. 2. The quasi-momentum < > vs. time for three 
different initial displacements: 40, 80 and 90 sites. When 
< k > reaches n/2 (i.e., for an initial displacement greater 
than ^cr ~ 84 calculated with Eq.(^) the system becomes 
modulationally unstable. 
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FIG. 3. The modulus square of the order parameter 
defined in Eq. (^) as a function of time, for three different 
initial displacements (40, 80 and 90 sites) and with the same 
parameters as in Fig. |^. When the quasi-momentum < k > 
reaches 7r/2 (i.e., for an initial displacement greater than ^cr), 
the order parameter drops to ~ 0; cf. Fig. ti. 
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FIG. 4. The center of mass vs. time for a supercritical 
initial displacement of 150 sites. Solid line: Gross-Pitaevskii 
equation; dashed line: discrete nonlinear Schrodinger equa- 
tion. The (displaced) center of the trap is at x = /im. 
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